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Abstract. We present CO(l-O) and CO(2-l) maps of the starburst/Seyfert 1 galaxy NGC 2782 obtained with the 
IRAM interferometer, at 2'.' 1x1'.' 5 and 0'.'7x0'.'6 resolution respectively. The CO emission is aligned along the 
stellar nuclear bar of radius ~1 kpc, configured in an elongated structure with two spiral arms at high pitch angle 
~ 90°. At the extremity of the nuclear bar, the CO changes direction to trace two more extended spiral features 
at a lower pitch angle. These are the beginning of two straight dust lanes, which are aligned parallel to an oval 
distortion, reminiscent of a primary bar, almost perpendicular to the nuclear one. The two embedded bars appear 
in Spitzer IRAC near-infrared images, and HST color images, although highly obscured by dust in the latter. 
We compute the torques exerted by the stellar bars on the gas, and find systematically negative average torques 
down to the resolution limit of the images, providing evidence of gas inflow tantalizingly close to the nucleus of 
NGC 2782. We propose a dynamical scenario based on numerical simulations to interpret coherently the radio, 
optical, and molecular gas features in the center of the galaxy. Star formation is occurring in a partial ring at 
~ 1.3 kpc radius corresponding to the Inner Lindblad Resonance (ILR) of the primary bar; this ring-like structure 
encircles the nuclear bar, and is studded with Ha emission. The gas traced by CO emission is driven inward by the 
gravity torques of the decoupled nuclear bar, since most of it is inside its corotation. N-body simulations, including 
gas dissipation, predict the secondary bar decoupling, the formation of the elongated ring at the ~lkpc-radius 
ILR of the primary bar, and the gas inflow to the ILR of the nuclear bar at a radius of ~200-300 pc. The presence 
of molecular gas inside the ILR of the primary bar, transported by a second nuclear bar, is a potential "smoking 
gun"; the gas there is certainly fueling the central starburst, and in a second step could fuel directly the AGN. 

Key words. Galaxies: individual (NGC 2782) - Galaxies:starburst - Galaxies: spiral - Galaxies: kinematics and 
dynamics — Galaxies: interstellar matter 



1. Introduction 
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* Based on observations carried out with the IRAM 
Plateau de Bure Interferometer. IRAM is supported by the 
INSU/CNRS (France), MPG (Germany), and IGN (Spain). 



Molecular gas is the dominant gas component in the inner 
regions of spiral galaxies, making CO lines unique tracers 
of nuclear gas dynamics. As such, they are also a pow- 
erful diagnostic for identifying how active nuclei (AGN) 
are fueled. To feed an AGN through accretion, there must 
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be an adequate supply of gas whose angular momentum 
has been sufficiently reduced to enable inflow within the 
small spatial scales surrounding the black hole. Although 
there is rarely a lack of circumnuclear fuel, it is not yet 
clear how angular momentum is removed to enable nuclear 
accretion. 

To better understand how AGN are fed and main- 
tained, we have been conducting for several years now 
a high-resolution h i gh-sen sitivity CO survey (NUGA, 



Garcfa-Burillo et al.l l2003h of galaxies at the IRAM 
Plateau de Bure Interferometer (PdBI). Altogether we 
have observed 12 galaxies in two CO transitions with 
up to four configurations of the array, so as to achieve 
the most sensitive (typically ~ 2 — 4 mJybeam"^ in 
lOkms"^ channels) and the highest resolution (1-2") sur- 
vey currently available. The results of the NUGA sur- 
vey so far have been surprising: there is no single un- 
ambiguous circumnuclear molecular gas feature connected 
with the nuclear activity. On e- and two-armed instabilities 
( Garcfa-Buri llo et al. ''2003'), well-ordered rings and nu- 
clear spirals (iCg mb es et al. 2004') , circumnucl ear asymme- 
tries ( Krips et al.ll2005f) . and large-scale bars ( Boone et al 



20071) are among the variety of molecular gas morpholo- 
gies revealed by our survey. Moreover, an analysis of 
the torques exerted by the stellar gravitational poten- 
tial on the molecular gas in four galaxies suggests that 
the gas tends to be driven away from the the AGN 
( ^ 5 pc), rather than toward it (|Garcfa-Burillo et al. 
2005h . Nevertheless, these dynamics do not correspond 



to the violent molecular outflows and superwinds pre 



dieted in AGN feedback inodels (e.g., iNaravanan et al 



2006t Hopkins fc HernquistI 2006[ ). because the observed 
velocities are much too small. 

Much of the explanation of this va riety of morphologies 



appe ars to be related to timescales (jGarcfa-Burillo et al 



20051) . It is well established that large -scale bars trans- 
port gas inward very effic iently (e.g.. Combes fc GerinI 
1985t Sakamoto et al.lT999f l. and there is ve ry little doubt 
that bars can drive po werful starbursts ( Knapen et al.l 



2002t iJogee et al.l l2005fl . However, no clear correlation 



betwe en bars and nuclea r acti vity has yet been found 
(e.g., iMulchaev fc ReganI 119971 ). This may be because 
the timescales for bar-induced gas inflow and AGN 
duty cycles are very different. Bars drive inflow over 
timescales ( ^ 300 Myr) that are similar to the typical 
gas-consumption timescales of a few times ^ 10^ yr found 
in nuclear starbursts fe.g.. iJogee et al.l 120051 ). But AGN 
accretio n-rate duty cycles are much shorter than this (~ 1- 



10 Mvr. iHeckman et al.lbooi iHopkins fc Hernquistll2006l : 



King fc Pringlell2007h . and there are several indications 
that active accretion occurs only intermit t ently over the 
lifetime of a galaxy dFerrarese et al.l 2001 : Marecki et al.l 



20031: IJaniuk et all I2OO4I : iHopkins fc HernquistI I2OO6I: 



King fc Pringlell2007f ). The resulting implication is that 



most AGN are probably between active accretion episodes. 



and catching galaxies with nuclear accretion "switched 
on" may be difficult. 

In this paper, the ninth of the NUGA series, we present 
observations that suggest that we have found one of these 
potentially rare AGN with possible gas inflow in the 
current epoch. NGC 2782 is an early-type spiral galaxy 
[SABa(rs)] with peculiar morphology. In addition to a pro- 
nounced stellar tail or she et ^20kpc t o the east, it shows 
three optical ripples (e.g.. Smith! 1994) thought to be sig - 
natures of tidal interactions (|Schweizer fc Seitzeii Il988[ ) . 
A massive Hi plume extends ~54kpc to the northwest, 
and the neutral atomic gas in the inner disk is coun- 
terrotati ng with resp ect to the gas motions in the outer 
regions (|Smithl [l99lh . The central regions of NGC 2782 
host a massive nuclear starburst, with a far-infrared 
(FIR) luminosity of 2 x 10^'^ Mq, comparable to that in 
M82 (Deyereux 1 989). Three-dimensional optical spec- 
troscopy ( Yoshida et al.lll9^9h shows evidence for a high- 
speed ionized gas outflow, with the bipolar structure in 
the radio continuum indicative of a confined superbubble 
(jjogee et al.l Il998l ) . In the outflow, there are also high- 
excitation extran uclear emission li nes thought to be due 
to shock heating (|Boer et al.lll992r ). 

Until recently, it was thought that the outflow and 
energetics in NGC 2782 were powered by a starburst 
alone, but recent radio and X-ray observations reveal 
an optically-hidden AGN. MERLIN and EVN/VLBI ob- 
servations show a high-brightness-temperature extremely 
compact (^0"05) radio so urce, unambiguous evidence of 
an AGN (|Krips et al.ll2007^ . NGC 2782 is also a Compton- 
thick X-ray source wit h a 6.4 keV iron feature coming from 
its innermost regions (jZhang et al.,.2006,) . There is a bright 
unresolved X-ray core and extended emission roughly co- 
incident with the radio morphology in the high-resolution 
(- 1") image by^aikia et al. (1994). 

The nuclear region of NGC 2 782 has been obse rved pre- 
viously in the 12 CO (1-0) line bv llshuzuki I (Il9 94') with the 
Nobeyama Millimeter Array and bv lJogee et al.. (1999) us- 
ing the Owens Valley Radio Observatory. We reobserved 
NGC 2782 in 1^00(1-0) at PdBI with better spatial res- 
olution and a sensitivity roughly five times that of previ- 
ous observations, and for the first time in the ^^C0(2-l) 
line. This enables a rigorous derivation of the torques act- 
ing on the molecular gas in NGC 2782, and a quantitative 
assessment of the infall of material to the nucleus. We 
first present our new observations in Section |2l together 
with our multiwavelength imaging dataset. The morphol- 
ogy and kinematics of the molecular gas are discussed in 
Section |3l and we describe the stellar structure in Section 
m and the starburst episode in Section |5l We then derive 
the gravitational potential from the infrared image, and 
infer the torques acting on the molecular gas in Section jS) 
Finally, we present numerical simulations which motivate 
our proposed scenario of decoupled double bars in this 
galaxy. The molecular gas in NGC 2782 is apparently be- 
ing driven inward by the nuclear bar, decoupled from the 
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primary bar, since we detect azimuthally averaged torques 
which are negative down to the resolution limit of our im- 
ages. 

2. Observations 

We observed NGC 2782 with the IRAM PdBI in the 
ABCD configuration of the array between December 2001 
and February 2003 in the 12CO(1-0) (115 GHz) and the 
i2CO(2-l) (230 GHz) rotational transitions. The PdBI re- 
ceiver characteristics, the observing procedures, and the 
image reconstruction are t he same as those described in 
Garcfa-Burillo et al.l (j2003l ). The quasar 3C273 was used 
for bandpass calibration and the quasar 0923-1-392 was 
used to calibrate both the phase and the atmospheric vari- 
ations. We used uniform weighting to generate 2-1 maps 
with a field of view of 21" and natural weighting to pro- 
duce the 1-0 maps with a field of view of 42". Such a 
procedure maximizes the flux recovered in CO (1-0) and 
optimizes the spatial resolution in C0(2-l). 

The rms noise a in lOkms^^ wide velocity chan- 
nels is 2.0mJybeam~^ and 5.2mJybeam~^, with beam 
sizes of 2"lxl"5 and 0"7x0"6 at 115 and 230GHz, 
respectively. At a level of ^ 3cr, no 3 mm (1mm) 
continuum is detected toward NGC 2782 to a level of 
ImJybeam"^ (3mJybeam~^). The conversion factors be- 
tween Jybeam^^ and K arc 30KJy^^beam at 115 GHz, 
and 58KJy~^beam at 230 GHz. By default, all veloc- 
ities are referred to the heliocentric recession velocity 
uo = 2545kms^^ and (Aa, A(5) offsets are relative to 
the phase tracking center of the observations (RA2000J 
Dec2ooo) = (09''14"05.08^ 40'*06™49.4''). The displayed 
maps are not corrected for primary beam attenuation. 

We will assume a distance to NGC 2782 of D =35 Mpc, 
wh ich is derived from t he local velocity field model given 
in Mould et al. ( 2000( ) and a Hubble constant Hq — 
73kms~^ Mpc~^. At this distance 1" corresponds to 
171 pc. 



2.1. Optical and infrared images 

We retrieved HST archival images of NGC 2782 with 
WFPC2 in the F555W and F814W filters. Cosmic rays 
were eliminated, and the images were calibrated a nd con- 
verted to V and / as described in lHoltzman et al] ([1995). 
We performed an astrometric calibration using stars from 
the U.S. Naval Observatory Astrometric Catalog Bl.O 
(USNO-Bl.O). Five stars from this catalogue appear in 
the 1600X 1600 WFPC2 image of NGC 2782 and they were 
used to derive the astrometric solution with imwcs in the 
WCSTools packageQ. The solution has an rms uncertainty 
of 0"24, or 2.4 WFPC2 (mosaic) pixels, and differs from 
the original HST one by ^ 1". We made V—I color images 



^ Available from 
\http://tdc-w'ww. harvard. edu/software/wcstools\ 



by converting the flux units to magnitudes and subtract- 
ing the two magnitude images. 

IRAC images at 3.6 /im were retrieved from the 
Spitzer archive. We started with the Basic Calibrated 
Data im ages, and aligned and co mbined them with 
MOPEX (jMakovoz fc Marleaul [20051 ) which accounts for 
distortion and rotates to a fiducial coordinate system. 
1"20 pixels were imposed for the final image, roughly the 
same as the original IRAC detector. Significant banding 
was present from a bright star in the field, and this was 
corrected for by interpolation over the affected rows be- 
fore combination with MOPEX. Combination of these cor- 
rected images showed no discernible effect of the original 
flaw. 

The infrared and optical images will be discussed in 
Sections H] and O 

3. Molecular gas results 

Figures [T] and [2] show the velocity-channel maps of 
^^CO(l-O) and "'^^00(2-1) emission in the central region 
of NGC 2782. The kinematics show the typical signature of 
a rotating disk, together with non-circular motions more 
prominent to the south. The western (eastern) side of the 
CO bar is red- (blue-) shifted relative to the systemic ve- 
locity Vsys, which is fltted to be at 2555±10kms~^. This 
agrees very well wit h the Hi heliocentric velocity of 2555 
kms-^ (|Smithlll994l ). 

The dynamical center was derived from the center of 
symmetry of the velocity fleld, and chosen to maximize 
the velocity dispersion. It is determined to lie at (02000 j 
(^2000)= (09''14'"05.11", 40°06'49.24"), not far from the 
phase-tracking center. The position of the dynamical 
center coincides alm ost exactly with th e compact radio 
source discovered bv iKrips et al. ( 2007[ ): (0:2000 j (52000) = 
(09''14™05.11^ 40°06'49.32"). Both are also within the er- 
rors of the position o f the "C" co mponent of the NGC 2782 
nucleus described bv lSaikia et al. (1994). Therefore we as- 
cribe the position of the AGN and the dynamical center 
of NGC 2782 to the position derived here from the CO 
kinematics. 



3.1. CO morphology 

The overall distribution of the molecular gas is illus- 
trated in Figure [3] which shows the velocity-integrated 
CO intensity maps, achieved by integrating channels from 
V = -230 to 230kms-i. The naturally weighted 1^00(1- 
0) map is shown in the top panel of Fig. [31 and the 
uniform-weight ^^C0(2-l) map in the bottom panel. After 
correcting for primary beam attenuation and considering 
only measurements with a signal-to-noise of 3cr or greater 
on the data cube, our naturally weighted (1-0) observa- 
tions in the central 42" detect an integrated CO flux of 
150 Jy kms~^. This corresponds to 65% of that in the Five 
College Radio Astronomy Observatory single dish mea- 
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surements (within 45": 230 Jvkms~^. lYoung et al.ll995 ) 
and r oughly 78% of the flux measured by Ijogee et al 
( I999I 195 Jykms"^). Our measurements are sensitive to 
small-scale structure but are missing some fraction of the 
diffuse compon ent. We use t he sh ort-spacing single dish 
observations bv lYoung etHI (|l995h to derive the gas mass 
for numerical simulations in Sect. [7731 




Aa(arcsec) 



Fig. 3. Top panel: a) The natural-weight map in ^^C0(1~ 
0) shown in grey scale with contours ranging from 0.5 
to 8.0 Jykms^^ in 1 Jykms"^ intervals. Bottom panel: 
b) The uniform-weight map in ^^C0(2-l) shown in grey 
scale with contours ranging from 1 to 7.0 Jykms"^ in 1 
Jykms^^ intervals. 

The CO (1-0) emission is distributed i n a nuclear 
elongated structure already identified by IJogee et"al] 



(|1999[ ). Our observations clearly delineate the diffuse spi- 
ral arms extending to the north and south. The spirals are 
not symmetric, and there is substantially more emission 
to the southeast than to the northwest. The filaments at 
^10" to the north and south in CO(1~0) correspond to 
the spiral arms from the outer stellar oval (which we dis- 
cuss in n iore detail b e low). Indeed, the CO "spur" labeled 
"02" bv lJogee et al. ( 19991 ) appears to be associated with 
the extended southern spiral arm shown in our maps (see 
Fig.©. 

The central structure in the ^"^€0(2-1) map is clearly 
resolved; the higher resolution afforded by our new ob- 
servations shows a clearly elongated structure with spiral 
arms commencing at the ends of the elongation. The gas 
in the inner spiral arms is aligned along the elongation 
with a pitch angle of ^ 90°, which makes the feature re- 
semble more a bar than an inclined or edge-on disk. At 
the end of the feature, the gas changes direction to follow 
the outer spiral arms, which are situated at a lower pitch 
angle. Moreover, the molecular gas is clearly responding 
to the stellar oval/bar which we will discuss in Section 21 
Hence, we will refer to the circumnuclear molecular struc- 
ture as the "nuclear gas bar" , as distinct from the nuclear 
stellar bar (see below). 

The molecular gas mass within the 42" PdB primary 
beam field is = 1.4 x 10^ Mq, assuming the CO-to- 

H2 conv ersion factor X = 2 2 x 1 0^° cm~^(Kkms~^)~^ 
given mon fc Barrett! (|l99ll ) . Including the helium 

mass in the clouds (multiplying xl.36) gives M^f^+He = 
1.9 X 10*^ Mq. Most of the molecular gas mass in NGC 2782 
is in the nuclear gas bar-|-spiral, making it an extremely 
massive structure. Indeed, the circumnuclear molecu- 
lar gas component in NGC 2782 is particularly massive, 
roughly 3 times more so than most of the NUGA galax- 
ies studied so far (NGC 4826: iGarci'a-Burillo et all 12003 



NGC 7217: ICombes et aP \2odi NGC 3718: iKriBselal 



2005 NGC 4579: Garcia-Burillo et al. 2005, NGC 6951 



Garcia-Burillo et al.l 2005 all have molecul ar masses on 
the o rder of ~ 3 x 10^ Mq). Only NGC 4569 (iBoone et al. 
2OO7I ). with = 1.1 X 10^ Mq, is roughly comparable 



with NGC 2782. 



3.2. Kinematics 



Figure m shows the mean velocity field derived from C0(1- 
0) (top panel) and CO (2-1) (bottom). The star marks the 
position of the AGN (coincident with the dynamical cen- 
ter). The kinematic signature of a rotating disk is clearly 
seen in Fig.lU and in the body of the circumnuclear molec- 
ular spiral, there are few non-circular motions. However, 
the southern spiral arm in the CO (1-0) map (top panel of 
Fig. [4]) shows clear streaming motions, and appears to be 
somewhat decoupled from the nuclear spiral because of the 
velocity discontinuity toward the southeast. In the bottom 
panel of Fig. (H systematic kinks appear in the CO (2-1) 
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Fig. 4. 12CO(1-0) (a: top panel) and 12C0(2-1) (b: bot- 
tom) isovelocities contoured over false-color velocity maps. 
The AGN (dynamical center) position is marked with a 
star. 

We derive a kinematic major axis of PA = (75 zt 5)°, 
consi stent with previous determinations ('75°. |Jogee et al. 
1999[l . We will assume that the inclination angle of the 
nucle ar spiral in NGC 2782 is ^30° (see also IJogee et"aL 
199i). A smaller inclination would result in sin i corrected 
velocities which would be too large to be consistent with 
even the largest spiral galaxies and ellipticals. A larger 
inclination seems equally unlikely because of the relatively 
straight extended spiral arms (see discussion in Sect.[S]); if 



the inclination were ^30°, the arms would be apparently 
more compressed relative to the line of nodes. That the 
galaxy disk is inclined at 30° is confirmed by the two- 
dimensional bulge-disk decomposition and the elliptical 
isophote fitting described in Sect. ID 
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Fig. 5. Top panel: a) Position-velocity diagram of 
i2CO(l-0) along the kinematic major axis of NGC 2782 
(PA = 75°) contoured over a grey-scale representation. 
Velocities have been rescaled to Fsy^ = 2555km s~^, and 
offsets are relative to the dynamical center. Bottom panel: 
b) The same for 12C0(2-1). 

Position- velocity (PV) diagrams along the major kine- 
matic axis of NGC 2782 are shown in Fig. O and along 
the minor axis in Fig. [6l In both figures, ^^CO(l-O) is 
given in the top panel, and ^^C0(2~I) in the bottom. 
The top panel of Fig. [5] reveals regular circular rota- 
tion together with mild streaming motions in the molec- 
ular gas. However, the rough alignment of the molecular 
feature at PA^88° with the line-of-nodes of the galaxy 
(PA~ 72 — 73°, see Sect. |4]) inhibits a clear signature of 
non-circular streaming motions. Consequently, the major- 
axis PV diagram cannot provide conclusive evidence for 



6 



Hunt et al.: NUGA: IX. NGC 2782 



300 



100 



A:l-0 




5) 


1 

cr 




1 


1 


1 



100 



B:2-l 


)7 : 


; - K 


<1> - 





Ay'(arcsec) 



Fig. 6. Top panel: a) Position-velocity diagram of 
^^CO(l-O) along the kinematic minor axis of NGC 2782 
(PA = 165°) contoured over a grey-scale representation. 
Bottom panel: b) The same for ^200(2-1). 



out-of-planc kinematics or for kinematic decoupling of the 
nuclear structure from the larger-scale stellar disk. The 
kinematics shown by the minor-axis PV in Fig. [5] are quite 
regular in the extended regions, but show significant ve- 
locity dispersion close to the nucleus. This could reflect 
an unresolved rotational velocity component, suggesting 
a large central dynamical mass (see below). Finally, Fig. [5] 
shows no evi dence for molecular gas outflow, as tentatively 
suggested bv lJogee et al. I ()l999l) . 



3.3. The rotation curve and dynamical mass 

We have derived a rotation curve (RC) from the PV dia- 
gram taken along the kinematic major axis of NGC 2782 
at a PA = 75 ° . The terminal velocities were derived by 
fitting multiple Gaussian profiles to the spectra across 
the major axis. The fitted velocity centroids, corrected 
for sini [i = 30°), give Kot for each galactocentric dis- 
tance. Although for both lines the RCs to the south (nega- 



tive velocities) are slightly steeper within the inner 300 pc 
than those to the north (positive), we averaged together 
the two curves derived from either side of the major axis. 
Hence, the resulting Vrot may be slightly shallower in the 
inner regions than the true mass distribution would imply. 
Because data for both lines were consistent, we combined 
both curves into an average by spline interpolation. The 
final RC is shown in Sect. l7.2l in the context of our numer- 
ical simulations. 

The peak velocity is ~ 170/ sin i kms^"'^, or 
^^340 km s^^ for i = 30°, obtained at a radius of 850 pc 
(5") to the south, but at ^1 kpc (~6") in the RC obtained 
from averaging both lines and both sides of the major axis. 
From the peak velocity, we can derive a rough estimate of 
the dynamical mass: M{R) = 2.32 x l(f>RV'^{R) where 
M{R) is in Mq, R in kpc, and V in kms^^. Assuming the 
most flattened disk-like distribution (i.e., including a mul- 
tiplicative constant of 0.6), at a radius of Ikpc, we would 
?mdMdyn = 1.6xlOi''M0. Within a 1.7 kpc radius, where 
virtually all of the observed molecular gas is located, as- 
suming a flat rotation curve (see Sect. 17. 2p . we would in- 
fer Mdyn = 2.7 X 1O^°M0. Comparing with the mass in 
molecular gas including helium would give a molecular 
mass fraction of ~ 7% in this region, in agreement with 
Jogee et all (|l999[ ). 



4. Stellar structure 

Figure[7]shows the CO(l-O) and (2-1) total intensity maps 
overlaid on the IRAC 3.6 /im image described in Sect. 12.11 
This wavelength traces very well the massive component 
of the stellar populations in galaxies, and has the added 
advantage of very low extinction even compared to the 
K band. Two main inner structures can be seen in Fig. 
[71 an inner stellar bar /oval with a PA^88° and diame- 
ter of '^15", and an outer oval at PA~10° and diame- 
ter of '^30'o Previous K-haad observations (IJogee et al 



19991 ) identified similar features, but with slightly different 



PA's and a larger radius for the outer oval than we find 
here (25" vs. ~15"). Both structures are clearly present 
in NGC 2782, and as we shall see in Sect.O contribute to 
the dynamical perturbations in this galaxy. 

A large-scale view of NGC 2782 at 3.6 /im is illustrated 
in (the left panel of) Fig. [H The detached stellar sheet or 
tail outside the main galaxy disk is evident to the east, 
and the elongation or distortion of the disk to the west 
i s clearly vis i ble. There are als o stellar arcs or "ripples" 
(|Smithlll994 Ijogee et al.lll999f) about 25" to the west of 
the nucleus. All these features are seen in optical image s 
(|Smithlll99i I Jogee et al.lll999f l and in Hi (|Smithlll99lh . 
Moreover, in Hi, there is a long (~ 54 kpc) plume or tidal 
tail extending to the north, with its origin at the western- 
most edge of the distorted disk ( Smith E99l[ ). 



^ These values have been derived by fitting elliptical 
isophotes to the image with the task ellipse in the 
IRAF/STSDAS package. 
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Fig. 7. Top panel: a) -'^^CO(l-O) total intensity contoured 
over a false-color representation of the IRAC 3.6 /im image 
of NGC 2782. The nuclear bar/oval and the outer oval can 
clearly be seen at PA^88° and '^10°, respectively. Bottom 
panel: b) The same for ^^C0(2-l). In both panels, the 
AGN (dynamical center) position is marked with a star. 



4.1. Bulge/disk decomposition 

To better understand the mass distribution and investi- 
gate the non-axisymmetric structure in the stars, we per- 
formed a two-dimensional bulge/disk decomposition on 
the IRAC 3.6 ^m image of NGC 2782 with qalfit, the pub - 
licly available algorithm developed bv IPeng et al. ( 2002 ). 
The MOPEX point-response function was used for convo- 
lution with the image, and the background sky level was 
fixed to the measured value rather than fit. Initial param- 
eters for bulge and disk were guessed by using the scaling 



relatio ns given bv iMoriondo et al" ( 1998f l and Hunt et al 
(j2004[ ). The final best fit was achieved with a generalized 
exponential (Sersic) bulge, an exponential disk, and a nu- 
clear point source. 

We ran several sets of fits, in order to experiment with 
masking the stellar tails or sheets predominant to the east. 
Fitting the whole image (~3.4'diameter) results in a bulge 
with shape index n = 3 (a de Vaucouleurs bulge has 
n — A), an effective radius of 1.1 kpc and an apparent 
inclination of 32°. This bulge contains about 60% of the 
total 3.6 fim light. The resulting disk, with a scalelength of 
2.9 kpc and an inclination of 41°, together with a nuclear 
point source contribute about 33% and 7%, respectively, 
to the total luminosity. The fitted aspect ratio of the disk 
is very close to that given for NGC 2782 in NED, which 
corresponds to a system inclination of 42° . The bulge PA 
of 19° is ill-determined because of its low apparent flat- 
tening, but the fitted orientation of the disk (PA = 73°) 
is consistent w i th the angle of the line-of-nodes estimated 
bv lJogee et aP (|l999h and with the kinematical major axis 
found in Sect l3.2l 

Masking the stellar "sheets", and confining the fit to 
the undisturbed portion of the outer regions (within a 
^^2. 4' diameter, 24 kpc) gives a slightly different, probably 
more reliable, fit. The bulge has a steeper shape index 
(n = 4), the nuclear point source is smaller, and the disk 
is less inclined. The fraction of bulge-to-total luminosity 
remains the same, ~60%. The masked fits give a disk incli- 
nation of 33°, more similar to that of the bulge, and more 
consistent with the visual aspect of the galaxy. Indeed, in- 
spection of Fig. [8] suggests a strikingly round appearance, 
at least in the regular portion of the disk out to a radius of 
J^50" (~9 kpc). This impression is confirmed by fitting the 
isophotes to ellipse^ which shows that at galactocentric 
distances as large as 12.5 kpc, the fitted system ellipticity 
implies an inclination of ^ 33° at PA = 72°. Hence, we 
confirm ^30° for the system inclination (as described in 
Sect. l3.2| ). obtained with the masked bulge-disk decompo- 
sition. 

The fit provides a convenient axisymmetric model for 
unsharp masking. The large-scale residuals from the fit de- 
scribed above are shown in the right panel of Fig. [8] The 
correspondence with the features in the unsharp masked 
optical images (iSmith 1994) is excellent. The stellar sheet 
to the east, and the ripples and distortion of the disk to 
the west, are clearly revealed. The small-scale residuals 
of the bulge/disk decomposition are shown in Fig. [9] with 
CO(l-O) (left panel) and C0(2-l) (right) overlaid in con- 
tours; both the "masked" fits and the fits to the entire im- 
age give virtually identical residuals in the circumnuclear 
region. The large (1'.'2) pixels of the IRAC image impede 
detailed comparison, but the residuals from the fit have an 
m = 2 structure, and strongly resemble a stellar bar. Some 
small component of residuals could perhaps be hot dust 



^ With the task ellipse in the IRAF/STSDAS package. 
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rather than stars, because of the red — L color excess 
obser ved in the center of NGC 2782 ( Hunt fc Giovanardil 
19921 ). Nevertheless, dust cannot contribute significantly 
to the structure of the residuals because the K — L color 
is redder outside the bar- like structure than within it. In 
fact, the position of the bar- like residuals coincides per- 
fectly with the CO emission, and provides a classic illus- 
tration of the theoretical behavior of the spiral response 
of a gas component to bar forcing. In this case, the gas is 
phase shifted in advance relative to the bar (leading), as 
can be seen to the northeast and the southwest in the "S" 
shape of the spiral. 



5. Star formation and dust in NGC 2782 

While stellar structure can be more readily inferred 
from near- infrared (NIR) wavelengths, dust extinction 
and star formation are best investigated in the optical. 
Figure \W\ shows the CO intensity maps overlaid on the 
HST/WFPC2 F814W image. The F814W emission of the 
circumnuclear region shows a strong asymmetric excess 
to the northwest of the AGN (hereafter called the "NW 
excess "), roughly coincid ent with the structure in the Ha 
maps ( Jogee et al. 19991 ). This excess is contiguous with 
an entire arc of bright knots to the north running from east 
to west, again reflecting the Ha emission-line morphology. 
Part of this arc lies within the ^■^ CO (2-1) emission, while 
the excess to the northwest is outside of it. There is also 
an extended F814W excess generally to the south, leading 
to a lopsided appearance at this wavelength. 

Figure [TT] presents the CO intensity maps overlaid on 
the HST/WFPC2 V-I image described in Sect. O The 
red V — I color delineating the spiral arm to the northwest 
is neatly traced by the CO(l-O) emission. This red arm 
connects to the CO spiral with a hook-like structure which 
winds around the NW excess, following the CO emission 
to the nucleus. The blue V — I colors of the F814W arc to 
the north of the AGN suggest that it is tracing the same 
star-formation episode revealed by Ha. 

This northe rn arc of star formation was identified 
by I Jogee et"aL 



and IJogee et al.l (|l999l ). and 



curs roughly cospatially with an outflow thought to be 
driven by the star burst. The outflow is s een at radio wave- 



lengths and in optical emission lines (IBoer et al 



Saikia et al. 1994|; iJogee et al. 1 119981 : lYoshida et al 



1992; 



1998; 



Jogee et al., 1999V It is relati vely north-south s ymmetric 
in the 6cm radio continuum ( Saikia et al.lll99^ . but has 
the form of a "m ushroom" in the o ptical (clearly illus- 
trated by Fig. 8 in iJogee et al.lll999[ ). The arc defines the 
"cap" of the mushroom, and the slightly narrower "stem" 
extends southward, ending at the blue knot roughly ~7" 
south and 1" east of the nucleus. This mushroom shape 
is not clearly seen in the V — I image. The roughly nor- 
mal - / colors (~0.8-l) to the south of the C0(2-l) 
emission suggest that the unreddened stellar population 
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Fig. 10. Top panel: a) 1^00(1-0) total intensity 
contoured over a false-color representation of the 
HST/WFPC2 F814W image of NGC 2782. Bottom panel: 
b) The same for 12C0(2-1). In both panels, the AGN (dy- 
namical center) position is marked with a star. 



of the bulge is obscuring the star formation to the south, 
a conclusion which we will motivate in the following. 

The three-dimensional orientation of the galaxy disk 
can be inferred on the basis of the red V — I colors in 
the extended spiral arm in the northwest. The roughly 
symmetric spiral arm to the southeast is not reddened, 
implying that the spiral structure to the northwest is be- 
ing seen through an obscuring dust lane, coplanar with 
the inclined stellar disk. If this is true, then the north 
would be the near side of the disk, and the south would 
be the far side. This orientation would also mean that the 
spiral arms are trailing, which is the most probable con- 
figuration. The implication is that, toward the south, the 
stellar disk is seen through the intervening bulge. The rel- 
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Fig. 11. Top panel: a) 12^0(1-0) total intensity 
contoured over a false-color representation of the 
HST/WFPC2 V -I image of NGC2782. The spiral arm 
to the northwest is clearly traced by the CO (1-0) emis- 
sion. Bottom panel: b) The same for ^^C0(2-l). In both 
panels, the AGN (dynamical center) position is marked 
with a star. 



atively small bulge effective scalelength of 5"— 6" means 
that ~50% of the bulge light lies within the region of the 
CO emission. To the south, the bulge stars with normal 
V — I colors are overwhelming the young blue stars in 
the star-formation event. Hence, the southern part of the 
"mushroom" structure is not observed in V — I, but in 
emission lines, since the bulge is not expected to contain 
a significant amount of dust, nor does it emit contami- 
nating line emission. To the south most of the bulge is in 



front of the intervening dust lane, while to the north, dust 
obscures more of the bulge because most of the bulge light 
is behind it. This interpretation would also help explain 
the F814W image with its apparent lopsidedness to the 
south. In any case, the dust must be clumpy, and inho- 
mogeneities in the dust distribution allow random regions 
behind the dust to shine through. 

The AGN (dynamical center) is located at the north- 
west border of a tiny red arc ^80 pc to the south- 
east ("parachute" -like, see Fig. [11]). In this region, there 
are broad Ha lines which show considerable north-south 
asymmetry, unlike the narrow comp onent of the emissio n 
lines which is relatively symmetric ( Yoshida et al. 1999f ). 



Interestingly, this red arc in the color map is at the same 
position as the velocity discontinuities in the C0(2-l) spi- 
der diagrams (see Sect. 13. 2p . Because of the feature's ex- 
tremely small size, we can not make any reliable state- 
ments, but speculate that it could be the color (and kine- 
matic) signature of a small-scale AGN outflow out of the 
plane of the disk. 

6. Gravitational torques on the molecular gas 

Streaming motions along the bar led lJogee et al. ( 1999f) to 



conclude that there could be gas inflow feeding the star- 
burst. However, streaming motions by themselves do not 
imply any transfer of angular momentum, nor can they be 
used to infer the direction of such transfer. Hence, we now 
examine the torques exerted by the stellar potential on the 
molecular gas in NGC 2782, in order to assess whether 
the gas is, on average, flowing in or out. Although our 
PdBI map recovers ~65% of that from single-dish mea- 
surements, in the following analysis what matters is the 
structured gas component. The diffuse large-scale molec- 
ular gas contributes equal quantities of positive and nega- 
tive torques on the gas orbits; hence it does not influence 
the net torque. 

There is no clear evidence for non-coplanarity of the 
CO nuclear bar-t-spiral and the large-scale stellar disk. 
The stellar bar is oriented at PA'^ 88°, and the kine- 
matic major axis at ^ 75°, somewhat hampering the iden- 
tification of non-circular motions. Nevertheless, were the 
two structures to be non-coplanar, we would expect to 
see some kinematic evidence which is not seen. Hence, 
in disagreement wit h the schematic model presented by 
Yoshida et al. ( 19991 ). in what follows, we assume that all 
the flattened structures (nuclear CO spiral, nuclear stellar 
bar, and outer stellar oval) lie within the same disk plane. 

To derive the torques exerted by the stars on the gas, 
we ad opt the method presented in ICarcia-Burillo et al.l 
([20051) which relies on a stellar potential $ derived from 
a NIR image and on the gas response as inferred from the 
CO maps. We assume that the total mass is determined 
by the stars, essenti ally invoking an extreme "maximum 
disk" solution (e.g., Kent! 1987 ). in which the gas self- 
gravity and the dark matter halo are neglected. We fur- 
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ther assume that the measured gas column density N{x, y) 
derived from the CO total intensity maps is a reliable es- 
timate of the probability of finding gas at (cc, y) in the 
current epoch. This is reasonable, namely that CO traces 
the total gas content, since atomic gas in the nuclei of 
galaxies is typically a very small fraction of the total gas 
mass. 

As emphasized in Garcfa-Burillo et al.l ( 2005[ ). a sin- 
gle key assumption underlies the validity of our estimate 
of angular momentum transfer in the gas: we implicitly 
assume that the gas response to the stellar potential is 
roughly stationary in the reference frame of the potential 
over a few rotation periods. This would not be the case 
in the presence of strong self-gravity in the gaseous com- 
ponents. However, otherwise, even in galaxies with sev- 
eral stellar pattern speeds at different radii in the disk, 
numerical simulations suggest that the gas response ad- 
justs its response to the dominant stellar confi guration at 
a given radius (see Garcia-Burillo et all l2005l and refer- 
ences therein). 

6.1. The stellar potential 

The stellar potential in NGC 2782 was derived from the 
IRAC 3.6 /im image by first rebinning to O'.'IS pixels, then 
deprojecting with a PA = 75° and an inclination of 30°. To 
account for the vertical mass distribution, the rcbinned de- 
projected image was convolved with an isothermal plane 
model of consta nt scale height^l/12th of the radial dis k 
scalelength (e.g.. lOuillen et al.lll994HButa fc Blockll200l[ ). 
Figure [T^] shows the deprojected IRAC image, with the 
^^CO(l-O) and ^^C0(2-l) total intensity maps superim- 
posed in contours. The orientations of the stellar outer 
oval and the nuclear oval/bar are clearly visible even in 
the deprojected image and are marked with a solid line. 

In the absence of extinction and stellar population 
gradients, the mass-to-light ratio can be assumed to be 
const ant across this deprojected (face -on) galaxy image 
(e.g.. lFreemanl[l995 Persic et al. 1996f ): the mass is there- 
fore directly proportional to the 3.6 /im light distribution. 
With this assumption, we derive the potential $ by solving 
the Poisson equation: — AttG p using a Fast Fourier 
Tran sform (FFT) te chnique on a 256x256 Cartesian grid 



each rn-Fourier amplitude by the normalized ratio Qm{R)'- 



Combes et al. I [1990 1. Beyond a radius of 3.4 kpc, 



as m 

is set to 0, thus suppressing any spurious m — A terms 
in the derived potential. This is adequate to compute <i> 
in the central regions of the galaxy sampled by the PdB 
CO(l-O) primary beam. 

Following ICombes fc Sanderd (|l98lL we then expand 
the potential in Fourier components: 



$(i?,6') = ^o{R) + ^^„i{R)cos[m{e - e,, 



(1) 



- m$,„/i?|Fo(i?)| 



(2) 



which is the amplitude of the m*'' harmonic of the force 
relative to the total axisymmetric component. These are 
essentially the harmonics of the angular derivatives of 
the potential, expressed as normalized Fourier ampli- 
tudes. Figure [13] gives the radial run of the highest-order 
Qm{R) coefficients. The nuclear (R~lkpc) and the pri- 
mary (R2±2.6kpc) bars are clearly evident as local max- 
ima (see left panel); moreover, the phase angle changes as 
would be expected for the ~ 90° transition from the nu- 
clear to the primary bar {Q2 shown as a solid line, right 
panel). 

6.2. Derivation of the torques 

The angular derivative of the potential ^{R,d) gives at 
each location the torque field per unit mass T{x,y) = 
xFy — yFx] by definition T{x,y) is independent of the 
gas distribution in the plane of the galaxy. The az- 
imuthal average of r, using N{x, y) as the weighting func- 
tion, corresponds to the global variation of the specific 
gas angular momentum transfer occurring at this radius 
{dLs{x,y)/dt\g). The 9 subscript indicates an azimuthal 
average, evaluated at a given radius R{x,y). For conve- 
nience of calculation on a pixelized image, we express the 
torques in Cartesian coordinates. 

The rotation direction of the gas in the galaxy disk 
determines the sign of T{x,y): positive (negative) if the 
torque drives the gas outward (inward) at {x,y). To cal- 
culate the effective variations of angular momentum in the 
galaxy plane, we then weight the torques by the gas col- 
umn densities N{x, y) derived from the ^^CO(l-O) and 2-1 
lines: t{x, y)N{x, y). Fig. ll4l shows the normalized version 
of these maps r' where we have divided t{x, y)N{x, y) by 
the azimuthal average of the gas density. The ' superscript 
indicates this normalized r, so that we can use dimension- 
less quantities (both in Cartesian coordinates and below 
in polar ones). Both panels in Fig. (TJ] show the typical 
"butterfly" diagrarn due to an m = 2 pertur bation (e.g. 
Buta fc Blo"c5l200ll: [Laurikainen fc Saloll2002[ ). 

The radial gas flow induced by the torques is estimated 
by azimuthally averaging the effective variations of angu- 
lar momentum density in the plane: 



r'(i?) 



JgN{x,y) X (xFy-yF,) 
IeN{x,y) 



(3) 



and use a polar Fourier transform method to isolate the 
various angular modes It is convenient to represent 



t'(R), by definition, represents the azimuthally-averaged 
time derivative of the specific angular momentum L^. 
Azimuthal averages are derived with a radial binning com- 
mensurate with the 1'.'2 resolution of the IRAC image. 
Similarly to the torque maps (e.g., Fig. [14]), the sign of 
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t'{R) defines whether the gas is subject to a net gain (+) 
or loss (— ) of angular momentum. Specifically, we assess 
the efficiency of the AGN fueling by estimating the aver- 
age fraction of the gas angular momentum transferred in 
one rotation period by the stellar potential. This is done as 
a function of radius by defining the dimensionlcss function 
AL/L: 



AL 



dL 



X Tr, 



T'iR) 



X Tr, 



(4) 



where T^ot is the rotation period, and Lg is assumed to 
be well represented by its azimuthal average, i.e., Lg = 
R X Vrot- The absolute value of L/ AL dictates how much 
time the stellar potential will need to transfer the entirety 
of the total gas angular momentum. Assuming that the 
gas response to the stellar potential is stationary with re- 
spect to the potential reference frame during a few rota- 
tion periods, a small value of AL/L implies that the stel- 
lar potential is inefficient at present. Figure [T5l shows the 
AL/L relative torques derived from CO(l-O) (top panel) 
and 2-1 (bottom). Because of the limited resolution in the 
IRAC image (1 pixel~200pc), we have masked the inner 
pixel. Beyond a radius of ^1.9 kpc, the ^^C0(2-l) emis- 
sion is less well sampled than CO (1-0) because of the 
smaller beam size ('^22"). The two CO lines are thus com- 
plementary, as CO (1-0) effectively traces the torques 
beyond this radius, but at a more moderate spatial reso- 
lution than CO (2-1). It is evident that both lines give 
similar results, in the sense that we see systematically neg- 
ative average torques in the inner few kpc. 

7. Understanding the torques in NGC 2782 

The ot her galaxies in the NUGA sample analyzed 
so far ( Garcia-Burillo ct all l2003l: ICombes et all 12004 : 
Garcfa-Burillo et al.,.200alBoone et al.ll2007t ) mainly show 
positive torques near the nucleus, with negative torques 
at varying amplitudes only for radii ^500 pc. Only in 
NGC 4579 are the average torques negative down to 
200 pc, but the spatial resolution of the torque map for 
that galaxy is sufficiently good that we can see that, in- 
side that radius, the gas is flowing away from the nucleus. 
NGC 2782 is the first galaxy in our sample that shows sys- 
tematically negative torques down to the resolution limit 
imposed by our torque image. Moreover, in the entire re- 
gion probed by our CO maps, nowhere is the gas subjected 
to average positive torques. The relatively poor resolution 
(1'.'2) of the potential map (IRAC 3.6 /um) makes it im- 
possible to assess the torque amplitude inside a radius of 
200 pc0. Indeed, the slightly positive torques in Fig. [15] are 
unreliable, and may simply be an artefact of insufhcient 



* As in lGarcfa-Burillo et al.l (|2005l ). we attempted to use the 
HST WFPC2 F814W image corrected for extinction. However, 
in NGC 2782, the extinction correction proved to be unreliable, 
and we were forced to rely only on the Spitzer/IRAC image. 
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Fig. 15. Top panel: a) The relative torque in dimension- 
less units of dL/L derived using the potential obtained 
from the Spitzer/IRAC 3.6 /im image (1.2" pixels, under- 
sampled beam). Bottom panel: b) The same for ^^C0(2- 
1). The ^^C0(2-l) emission is less well sampled than 
^^CO(l-O) beyond a radius of ^1.9 kpc, because of the 
smaller beam size (FWHP~22"). Despite the relatively 
poor resolution of the gravitational potential image, both 
CO transitions give similar results; we see systematically 
negative average torques within the entire region probed by 
our observations. 



spatial resolution. Nevertheless, the systematically nega- 
tive trend of the average torques in NGC 2782 hints that 
gas is being transported tantalizingly close to the AGN. 
This gas is certainly fuelling the nuclear starburst, and 
may, subsequently, be feeding the AGN. In the following, 
we develop a dynamical scenario which will explain why 
we think this occurs in NGC 2782. 

7.1. Main features 

Although the outer parts of NGC 2782 show clear signs 
of a past interaction (tidal tails, distorted disk, stellar 
sheets), the inner parts appear dynamically relaxed, not 
unexpectedly since the dynamical timescale is shorter by 
one or two orders of magnitude at r~ Ikpc, than at r~ 
20kpc. 

The previous sections have shown that NGC 2782 con- 
tains many features of a double-barred galaxy. Although 
the nuclear bar is unfortunately almost aligned with the 
galaxy line of nodes, the disk is nearly face-on (30°), which 
helps to distinguish the morphological features. In the 
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near-infrared, there is a nuclear stellar bar of radius ~ 
1 kpc, embedded within an outer stellar oval of radius ~ 
2.5 kpc. This outer oval, which we will call the primary 
bar, is oriented at PA ^10°, misaligned (and not quite 
perpendicular) with the nuclear stellar bar at PA '^88°. 
The molecular gas distribution is mainly elongated along 
the nuclear bar, and reveals two spiral features emerging 
from its ends, aligned with the conspicu ous dust lanes al- 
ready noted bv lJogee et al. (|l998l . ll999l) . These dust lanes 
are strikingly straight, reminiscent of barred galaxies with 
dust lanes parallel to the bar, but located on the bar's 
leading edge on either side of the nucleus. We therefore in- 
terpret the large-scale CO(l-O) features corresponding to 
the straight dust lanes as the gas response to the primary 
bar; they have indeed a similar orientation. The inner elon- 
gated molecular morphology, evident in both CO (1-0) and 
CO (2-1), constitutes the gas response to the nuclear bar. 

The Ha morphology is more puzzling, because of a 
cl umpy dust obscurat ion which masks part of it. As noted 
bv lJogee et al. ( 1998I) . there are two main features in Ha: 
one is the outflow driven by the starburst, perpendicu- 
lar to the nuclear st ellar bar, and alig ned with the radio 
continuum outflow (jSaikia et al.l 119941 ) . and the other is 
the elongated ring, of about 1 kpc in radius, encircling the 
nuclear bar. The ring is particularly demarcated on the 
northern side, where it forms a 180° arc, making the struc- 
ture resemble more a half-ring or what we called the "cap 
of the mushroom" in Sect. [5l But there is also a south- 
ern equivalent (dubbed the "mushroom stem"), largely 
obscured by dust. Toward the west, there is a remarkable 
dust finger, devouring the ring, which superposes quite 
well with the CO distribution (see Fig. fTTj). 

The two different Ha components are well distin- 
guished by their different kinematics: large outflowing ve- 
locities in the north-south component (bubble to the south 
in front of the dust lane, and obscured flow in the north be- 
hind it), and a no rmally rotating ring around the nuclear 
bar (jYoshida et a l. 1999). We interpret the low- velocity 
rotating component as the nuclear ring, associated with 
the inner Lindblad resonance (ILR) of the primary bar, 
and encircling the nuclear bar, i.e. close to the corotation 
(CR) of the nuclear bar (more exactly at the UHR, or 
4:1 resonance). The dust features are more conspicuous 
toward the north because the near side is the north side. 
The bulge is more evident toward the south, masking the 
dust obscuration from the disk's far side. But this does 
not imply that there is no dust and gas there, obscuring 
the southern half-ring in Ha: on the contrary, there is CO 
emission in this southern part, coinciding with holes in the 
Ha distribution. 

The Ha emission is the most intense in the partial 
ring, which lends support to our ILR interpretation, since 
it is well-known that resonances are locations where the 
gas is transiently stalled due to the torques provided by 
the stellar bar. Thus, the gas has time to form young 
stars, as typically revealed in many barred galaxies, such 



as NGC 2997, NGC 4314, or NGC 4321 (iMaoz et al.l ll996l: 
Benedict et"al1 l2002t lAUard et all 120061 ). In the case of 
NGC 2782, because of the secondary bar, gas is then driven 
further inward, as suggested by the CO straight spiral 
arms (pitch angle close to 90°) aligned with the nuclear 
bar. This gas also forms stars, explaining the presence of 
Ha also in the nuclear disk. 

7.2. Dynamical interpretation 

The dynamical scenario we propose to explain these main 
features is the following: first a primary bar forms, most 
likely triggered by the interaction that gave rise to the dis- 
turbed morphology in the outer parts of NGC 2782. The 
gravity torques from this primary bar drive the gas in- 
ward, and create a gaseous nuclear ring at its ILR. Intense 
star formation gives rise to the Ha nuclear ring, seen now 
as a sharp half-ring arc. In the me antime, the gas infiow 
has weak ened the primary b ar (e.g., Bournaud fc Combed 
[2002; Bo urnaud et al. 2005 ). and a nuclear bar has dy- 
namically decoupled, rotating at a different pattern speed. 
The two different pattern speeds are usually such that the 
corotation of the nuclear bar corresponds to the ILR of the 
primary bar, thus avoidi ng excessive chaos in the orbits 
( Friedli fc: Martinet 1993 ) . The nuclear bar then governs 
the gas infiow toward the center. Instead of being stalled 
at the primary ILR, now the gas is subject to the nega- 
tive gravity torques of the secondary bar, between its CR 
and its own ILR. This explains why we see CO emission 
inside the ILR of the primary bar, and why the torques 
computation reveals that gas is flowing inward toward the 
center. 

This scenario is slightly different from one where the 
gas decouples fro m the stars to form a gaseo us nuclear 
bar with a faster (Englmaier fc Shlosmanll2004l ) or slower 
(jHeller et al.ll200lh pattern speed than the primary bar. 
We propose here that the nuclear bar exists both in the 
stellar and gaseous components, and always with a higher 
pattern speed than the primary. 

7.3. Numerical techniques 

To better understand the physical mechanisms, and re- 
produce the observed morphology, we performed N-body 
simulations with stars and gas, including star formation. 
We believe that the nuclear bar decoupling is a disk phe- 
nomenon, occurring on short timescales, with little con- 
tribution of a slowly rotating spherical component, such 
as a bulge or dark matter halo, so we adopt the simplifi- 
cation of 2D simulations, with the spherical components 
added analytically. Self-gravity is included only for the 
disk (gas-l-stars) . 

The N-body simulations were carried out using the 
FFT algorithm to solve the Poisson equation, with two 
kinds of grids: Cartesian and polar, the latter optimiz- 
ing the spatial resolution toward the center. The useful 
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Cartesian grid was 256x256 (or 512x512 in total to avoid 
periodic images), with a total radius of 8 kpc. The cell size 
is then 62 pc. The softening was taken as 250 pc, in order 
to roughly reproduce the actual effective softening from 
the thickness of the disk. The polar grid is composed of 
NR= 80 radial, and NT= 96 azimuthal separations. The 
radial spacing of the grid is exponential at large scales, 
and linear at small scales, following the prescription by 
Pfenniger fc Friedlil (jl993l ) . The polar-grid cell size ranges 
from 10 pc at the center to ~lkpc at the outskirts (at 
22 kpc). The minimum softening is also fixed at 250pc, 
and is larger in the outer regions. The two different kinds 
of simulations, polar and Cartesian, were performed sep- 
arately, in order to gauge the effect of different spatial 
resolution on the results; both methods give similar out- 
comes. 

The stellar component is represented by 120000 parti- 
cles, and the gas component by 40000. Two rigid spherical 
potentials are added analytically, with a Plummer shape: 



$6Ji(r) = - 



' b.h 



where Mf, and rf, are the mass and characteristic radius 
of the stellar spherical component (bulge), and Mh and 
Th are the equivalent quantities for the dark-matter halo. 
These parameters are selected to fit the rotation curve of 
NGC 2782, and are reported in Table [H 

The stellar disk is initially a Kuzmin-Toomre disk of 
surface density: 

with mass M^;, rd = 3.25 kpc, and truncated at 7 kpc. It 
is initially quite cold, with a Toomre Q parameter of 1.2. 
The gas distribution is initially similar to the stellar one, 
with Qgas =1- The time step is 0.1 Myr. The initial con- 
ditions of the runs described here are given in Table [T] 
The component subscripts refer to the bulge (" 6" ) , and 
the halo (" W ) which are rigid, and the disk (" (i" ) which 
is "live". 



Table 1. Initial conditions parameters 



Run 




Mi, 


Md 




r gas 




kpc 


Mo 


Mo 


Mo 


% 


Run A 


0.2 


2.3el0 


1.4ell 


4.5el0 


10 


Run B 


0.5 


2.3el0 


g.OelO 


g.OelO 


10 


Run C 


0.2 


2.3el0 


1.4ell 


4.5el0 


5 


Run D 


0.5 


2.3el0 


g.OelO 


g.OelO 


5 



All masses are given inside 7 kpc radius; Vd and Vgas axe fixed 
at 3.25 kpc, and r^ at 7.5 kpc. 

The gas fraction Fgas is the fraction of the total mass (halo 
included) . 



The gas is treated as a self-gravitating component in 
the N-body simulation, an d its dissipation is treate d by a 
sticky particle code, as in ICombes &: Gerin (|l985l ). The 
initial gas-to-total mass ratio (Fgas) ranges between 5 
and 10%. Star formation is included, with a conventional 
Schmidt-law recipe, although this will not be discussed 
here, since the timescales are short enough that the gas de- 
pletion is not significant. At the present time, the observed 
gas mass in NGC 2782 within a ~3 kpc radius is about 
6 X 10^ MeE. The gas clouds are subject to inelastic col- 
lisions, with a collision cell size of 67 pc for the Cartesian 
code, and 117pc for the polar code (region where parti- 
cles are selected to possibly collide). This corresponds to 
a lower limit for the average mean free path of clouds be- 
tween two collisions. The collisions are considered every 
2 Myr. In a collision, the sign of the relative cloud ve- 
locities is reversed and the absolute values are reduced: 
relative velocities after the collision are only 0.65 times 
their original value. The dissipation rate is controlled by 
this factor. All gas particles have the same mass. 

The rotation curve corresponding to Run A is plot- 
ted in comparison to the empirical RC (see Sect. 13. 3p in 
Fig.mi All other runs have similar rotation curves. In the 
figure, we show the circular velocity Vdrc deduced from the 
potential (d$/dr = V^^^^/r). In the range of 200-800 pc, 
there are non-circular streaming motions (see Sect. 13.21) 
which contribute to the derived empirical VJ-ot- Indeed, 
we have calculated K-ot of the gas from the model, az- 
imuthally averaged at a given radius, and find streaming 
motions of amplitude similar to the observed discrepan- 
cies between the empirical RC as traced by CO and the 
predicted K:irc0. Given the uncertainties, and the asym- 
metry of the two sides of the major axis (see Sect. 13. 3t . all 
theoretical curves are compatible with the empirical ones. 
The data points for the epicyclic frequency resonances, 
marked by x and -I- in Fig. I16[ are calculated by numer- 
ically differentiating the smoothed empirical RC. These 
data are also generally compatible with the predictions of 
the numerical simulations. 

7.4. Simulation results 

Several runs have been carried out in order to test the 
various parameters, and evaluate the dependence of the 
results on the central mass concentration and on the ini- 
tial gas mass fraction. These parameters are crucial for 
the nuclear bar decoupling. Only four of these models are 
displayed in Table [T] 

We found that a high central mass concentration was 
ne cessary to decouple t he sec ondary bar, as already shown 
bv lFriedh fc Martinetl (|l993l ). This can be understood be- 



^ As derived f rom t he short-spacing single-dish observations 
of I Young et al] ijiggsh . see Sect. [3A] 

^ As noted in Sect. 13.31 the unfortunate coincidence of the 
PA of the bar and the disk line-of-nodes makes streaming sig- 
natures difficult to diagnose in PV and isovelocity diagrams. 
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Fig. 16. Rotation and derived frequency curves ft, VL—k/2 
and + k/2 adopted for the simulations (Run A), to- 
gether with the empirical curves derived from the C0(2- 
1) and (1-0) PV diagrams. The (green) open circles trace 
the empirical VJ-ot, coinciding with theoretical circular ve- 
locity Kirc (traced by a solid line) at R = 1 kpc and at 
small radii. The remaining solid lines are the empirical 
(black) and theoretical (green) curves. The long-dashed 
line is the Q. — n/A curve, and the short-dashed ones are the 
J7±k/2 ones. The corresponding data points are shown as 
X and the data have been masked within 500 pc to min- 
imize confusion. The two pattern speeds inferred from the 
simulations are shown as horizontal dotted lines, and reso- 
nances are marked with filled circles. The vertical dashed 
lines correspond to the nuclear bar length (~lkpc) and 
the ILR of the outer oval (^^l.Skpc), also roughly equal 
to the nuclear bar's CR at '^l.lkpc. The apparent mis- 
match of the empirical and theoretical RCs, in the range 
of 200— 800 pc, is due to non-circular streaming motions 
(Sect. 13111) and asymmetry in the two sides of the rotation 
curve (Sect. 13. 3p which contribute to the inferred average 



cause the two pattern speeds are close to the precessing 
rates of the supporting orbits, which are traced by the 
fl — k/2 curve. A high mass concentration produces a 
marked maximum of this curve in the center, with a large 
radial gradient favorable to the decoupling. Also, a high 
gas mass fraction of at least 5% of the total mass is re- 
quired for the decoupling. The nuclear bar is even stronger 
with a gas fraction of 10%. In Run B and D, with a scale- 
length of the bulge twice that of Run A and C, the mass 
concentration was not enough to produce the decoupling 
of the nuclear bar. 

Run A, with a small bulge scalelength, and a high gas 
mass fraction, provides the best fit for NGC 2782. The 
stellar and gaseous morphologies are displayed in Fig. [T7| 
at three epochs, shortly after the second bar decoupling. 



As seen in Fig. \T7l between ^700 and 750 Myr the gas is 
inflowing. Moreover, the gas is always "leading" the stars, 
and thus is slowed down by them. 

The Fourier analysis of the gravitational potential dur- 
ing the decoupling of the two bars clearly shows the devel- 
opment of the two modes. The Fourier analysis has been 
performed as a function of radius, every 2 Myr, during 
100 Myr in this period. The Fourier transform of the re- 
sults over the time axis gives the pattern speed as a func- 
tion of radius, as shown in Fig |18i with a resolution in Qp 
of 30kms~^ kpc^"'^. Of necessity, the temporal averaging 
must be performed over the lifetime of the transient fea- 
tures; this averaging "thickens" the value of Qp in the plot 
because of the time variation of the pattern speed of the 
nuclear bar. 

Figure [IH] lets us unambiguously isolate two bars in the 
potential, dominated by the stellar component: 



(a) The nuclear bar with fJf, ~270kms^^ kpc^^. Its coro- 
tation is at about 1.1 kpc, at a radius just slightly 
larger than its size of 1 kpc. 

(b) The primary bar with fit, ~ 65kms^^kpc^^ and an 
ILR at '^1.3 kpc, roughly coincident with corotation 
of the nuclear bar; the primary bar CR is at ^5 kpc, 
and its UHR at '^3.6-3.8 kpc. In the gas morphology, 
there is a (pseudo) inner ring coinciding with the UHR 
resonance (see Fig. [T7)) . 



From these identifications, we can interpret the elon- 
gated nuclear ring as the ILR of the primary bar at R ~ 
Ikpc, coinciding with the UHR/CR of the nuclear bar. 
This ring corresponds to the Ha arc, or half-ring, encir- 
cling the NIR nuclear stellar bar. Inside the ILR ring (of 
the primary bar), the gas morphology reveals two spiral 
arms with large pitch angle, now aligned with the nuclear 
bar. These gas arms are leading with respect to the stellar 
bar, and this phase shift means that the gravity torques 
from the stars are negative, driving the gas inward. These 
arms are winding up in a ring around the center, which 
corresponds to the ILR of the nuclear bar. Indeed, it can 
be seen from the frequency curves in Fig [16] that the pat- 
tern speed of the secondary bar (270kms~^ kpc~^) crosses 
the fl — k/2 curve at about 300 pc0. The two arms aligned 
along the nuclear bar can be identified with the CO fea- 
tures, observed with the interferometer. There is also a 
hole in the center in the gas distribution, which could cor- 
respond to the ILR of the nuclear bar. 



^ Indeed, the high maximum of the f2 — k/2 curve is cer- 
tainly overestimated by the epicyclic approximation, since the 
orbits are very elongated. Hence, although in our simulations, 
Qp of the secondary bar would indicate two ILRs, it is more 
likely that there is only a single ILR, roughly in the middle at 
~200pc. 
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7.5. Decoupled bars: gas infall in NGC 2782 

The previous simulations, adapted to the rotation curve of 
NGC 2782, support the scenario proposed in section mi A 
primary bar, perhaps triggered by the past interaction, has 
driven gas toward its ILR, and while this bar was weaken- 
ing, a secondary embedded bar has decoupled, prolonging 
the gas inflow inside the ILR ring, which is now conspicu- 
ous in star formation and Ha emission. The gas has been 
driven into two spiral arms, aligned with the nuclear bar, 
and is now winding up to the ILR of the nuclear bar, i.e. 
at about 200-300 pc radius. This scenario allows a coher- 
ent interpretation of the low-velocity Ha component and 
the entire CO morphology, as well as of our finding of 
negative gravity torques. In NGC 2782, there seems to be 
a "smoking gun", namely the presence of molecular gas 
inside the ILR of the primary bar, made possible by the 
decoupling of a secondary nuclear bar from the primary 
one. The gas there is certainly fueling the central star- 
burst, at the origin of the ionized gas bubble, and in a 
second step should fuel directly the AGN. The timescales 
in this central region are quite short, about 6Myr for the 
gas rotation period at 300 pc. Dynamical friction and vis- 
cous torqu es can then relatively quickly pro vide gas to fuel 
the AGN (jCarcia-Burillo et ahllBosl l2007t ). although the 
present simulations do not have enough spatial resolution 
to simulate that realistically. Alternatively, the ILR of the 
nuclear bar suggested by our models could impede further 
gas inflow inside the resolution limit of our images. 

According to our simulations, there are two salient 
characteristics of NGC 2782 which lead to the secondary 
bar decoupling, and the ensuing negative torques. The 
first is its conspicuous stellar bulge, as suggested by the 
bulge-disk decomposition described in Sect. 14.11 This is 
the biggest contributor to the central mass concentration, 
and the likely cause of the resonances that give rise to 
the gas buildup in the rings, with the subsequent decou- 
pling of the secondary nuclear bar. The second feature is 
the particularly high gas mass fraction in NGC 2782. As 
discussed in Sect. 13.11 the mass of the molecular gas in 
NGC 2782 is the most extreme of all the NUGA galaxies 
observed so far; within a region ~4kpc in diameter, the 
molecular gas comprises 7% of Mdyn- This centrally con- 
centrated molecular gas could be the result of the inter- 
action which NGC 2782 has undergone. A high gas mass 
fraction is a necessary condition for the decoupling of the 
secondary bar, and thus promotes infall. Finally, we could 
be observing NGC 2782 in a particularly favorable epoch 
for the fueling of its AGN. Future work on our statistical 
sample will help understand duty cycles and timescales 
for the creation and maintenance of active accretion in 
galactic nuclei. 



8. Summary and conclusions 

We summarize our main results as follows: 



(a) High-resolution ^^CO observations of NGC 2782 ob- 
tained with the PdB interferometer show emission 
aligned with the stellar nuclear bar of radius ~1 kpc. 
In CO (2-1), we clearly resolve two spiral arms at 
high pitch angle along the nuclear bar; in the ^^C0(1- 
0) maps, the gas changes direction to form the begin- 
ning of two straight dust lanes, aligned with an outer 
stellar oval reminiscent of a primary bar. 

(b) From the stellar gravitational potential inferred from 
the IRAC 3.6 /zm image and the gas density from 
the ^^CO distributions, we compute the azimuthally- 
averaged torques exerted by the stellar bars on the 
gas. This calculation reveals systematically negative 
torques in NGC 2782, down to the resolution limit of 
our image. 

(c) The negative torques and the dynamics inferred from 
our observations are explained by a scenario which is 
consistent with existing radio, optical, and molecular- 
gas features in the circumnuclear region of the galaxy. 
Star formation occurs in an Ha ring-like structure en- 
circling the nuclear bar at the ILR of the primary 
bar/oval. The gas traced by CO emission is infalling 
to the center as a result of gravity torques from the 
decoupled nuclear bar. The two high-pitch angle CO 
spiral arms are winding up toward the center in a nu- 
clear ring, corresponding to the ILR of the nuclear 
bar. 

(d) This scenario is supported by N-body simulations 
which include gas dissipation, and reproduce the sec- 
ondary bar decoupling, the formation of the elongated 
ring at the 1 kpc-radius ILR of the primary bar, and 
the gas inflow to the 200-300 pc-radius ILR of the nu- 
clear bar. 

The average negative torques revealed by our analy- 
sis show that infall to the central regions in NGC 2782 
is occurring in the present epoch. Our numerical simula- 
tions suggest that the negative torques producing inflow 
are caused by the decoupling of the primary and nuclear 
bar. This is made possible by the high central mass con- 
centration and large gas mass fraction in the galaxy, per- 
haps a result of the prior interaction which caused the 
stellar sheets, ripples, and Hi tidal tails in its outer re- 
gions. Future work will help understand what conditions 
are necessary for fueling an AGN, and how the gas dis- 
tribution, its kinematics, and the stellar gravitational po- 
tential interact to set the stage for gas inflow. 

Acknowledgements. LKH sincerely thanks the LERMA- 
Observatoire de Paris for generous funding and kind hospi- 
tality during the writing of this paper. We would like to 
thank the anonymous referee for insightful questions which 
helped clarify the text. This research has made use of the 
NASA/IPAC Extragalactic Database (NED) which is oper- 
ated by the Jet Propulsion Laboratory, California Institute of 
Technology, under contract with the National Aeronautics and 
Space Administration. 



16 



Hunt et al.: NUGA: IX. NGC 2782 



References 

AUard, E. L., Knapen, J. H., Peletier, R. F., & Sarzi, M. 

2006, MNRAS, 371, 1087 
Benedict, G. F., Howell, D. A., J0rgensen, I., Kenney, 

J. D. P., & Smith, B. J. 2002, AJ, 123, 1411 
Boone, F., Baker, A. J., Schinnerer, E., et al. 2007, A&A, 

471, 113 (NUGA VII) 
Boer, B., Schulz, H., & Keel, W. C. 1992, A&A, 260, 67 
Bournaud, F., & Combes, F. 2002, A&A, 392, 83 
Bournaud, F., Combes, F., & Semelin, B. 2005, MNRAS, 

364, L18 

Buta, R., & Block, D. L. 2001, ApJ, 550, 243 

Combes, F., Garcfa-Burillo, S., Boone, F., et al. 2004, 

A&A, 414, 857 (NUGA II) 
Combes, F., Debbasch, F., Friedli, D., & Pfenniger, D. 

1990, A&A, 233, 82 
Combes, F., & Gerin, M. 1985, A&A, 150, 327 
Combes, F., & Sanders, R. H. 1981, A&A, 96, 164 
de Jong, R. S. 1996, A&A, 313, 377 
Devereux, N. A. 1989, ApJ, 346, 126 
Englmaier, P., & Shlosman, I. 2004, ApJ, 617, L115 
Ferrarese, L., Pogge, R. W., Peterson, B. M., Merritt, D., 

Wandel, A., & Joseph, C. L. 2001, ApJ, 555, L79 
Freeman, K. C. 1992, Physics of Nearby Galaxies: Nature 

or Nurture?, 201 
Friedh, D., & Martinet, L. 1993, A&A, 277, 27 
Garcia-Burillo, S., Combes, F., Hunt, L. K., et al. 2003, 

A&A, 407, 485 (NUGA I) 
Garcfa-Burillo, S., Combes, F., Schinnerer, E., Boone, F., 

& Hunt, L. K. 2005, A&A, 441, 1011 (NUGA IV) 
Garcia-Burillo, S., Combes, F., Usero, A., & Gracia- 

Carpio, J. 2007, New Astronomy Review, 51, 160 
Heller, C, Shlosman, I., & Englmaier, P. 2001, ApJ, 553, 

661 

Janiuk, A., Czerny, B., Siemiginowska, A., & Szczerba, R. 

2004, ApJ, 602, 595 
Heckman, T. M., Kauffmann, G., Brinchmann, J., 

Chariot, S., Tremonti, C, & White, S. D. M. 2004, ApJ, 

613, 109 

Holtzman, J. A., Burrows, C. J., Casertano, S., Hester, 
J. J., Trauger, J. T., Watson, A. M., & Worthey, G. 
1995, PASP, 107, 1065 
Hopkins, P. F., & Hernquist, L. 2006, ApJS, 166, 1 
Hunt, L. K., & Giovanardi, C. 1992, AJ, 104, 1018 
Hunt, L. K., & Malkan, M. A. 1999, ApJ, 516, 660 
Hunt, L. K., Pierini, D., & Giovanardi, C. 2004, A&A, 
414, 905 

Ishuzuki, S. 1994, in lAU CoUoq. 140, Astronomy with 

Millimeter and Submillimeter Wave Interferometry, ed. 

M. Ishugiro & J. M. Welch (ASP Conf. Ser. 59; San 

Francisco: ASP), 292 
Jogee, S., Kenney, J. D. P., & Smith, B. J. 1998, ApJ, 494, 

L185 

Jogee, S., Kenney, J. D. P., & Smith, B. J. 1999, ApJ, 526, 
665 



Jogee, S., ScoviUe, N., & Kenney, J. D. P. 2005, ApJ, 630, 
837 

Kent, S. M. 1987, AJ, 93, 816 

King, A. R., & Pringle, J. E. 2007, MNRAS, in press 

(astro-ph/0701679) 
Knapen, J. H., Perez-Ramirez, D., & Laine, S. 2002, 

MNRAS, 337, 808 
Krips, M., Eckart, A., Neri, R., et al. 2005, A&A, 442, 479 

(NUGA HI) 

Krips, M., Eckart, A., Krichbaum, T. P., et al. 2007, A&A, 

464, 553 (NUGA V) 
Laurikainen, E., & Salo, H. 2002, MNRAS, 337, 1118 
Lequeux, J. 1983, A&A, 125, 394 

Makovoz, D., & Marleau, F. R. 2005, PASP, 117, 1113 
Maoz, D., Barth, A. J., Sternberg, A., et al., AJ, 111, 2248 
Marecki, A., Spencer, R. E., & Kunert, M. 2003, 

Publications of the Astronomical Society of Australia, 

20, 46 

Moriondo, G., Giovanardi, C, & Hunt, L. K. 1998, A&AS, 
130, 81 

Mould, J. R., Huchra, John P., Freedman, Wendy L., et 

al. 2000, ApJ, 529, 786 
Mulchaey, J. S., & Regan, M. W. 1997, ApJ, 482, L135 
Narayanan, D., Cox, T. J., Robertson, B., et al. 2006, ApJ, 

642, L107 

Peng, C. Y., Ho, L. C, Impey, C. D., & Rix, H.-W. 2002, 
AJ, 124, 266 

Persic, M., Salucci, P., & Stel, F. 1996, MNRAS, 281, 27 
Pfenniger, D., & FriedH, D. 1993, A&A, 270, 561 
Quillen, A. C, Frogel, J. A., & Gonzalez, R. A. 1994, ApJ, 
437, 162 

Saikia, D. J., Pedlar, A., Unger, S. W., & Axon, D. J. 

1994, MNRAS, 270, 46 
Sakamoto, K., Okumura, S. K., Ishizuki, S., & Scoville, 

N. Z. 1999, ApJ, 525, 691 
Schweizer, F., & Seitzer, P. 1988, ApJ, 328, 88 
Smith, B. J. 1991, ApJ, 378, 39 
Smith, B. J. 1994, AJ, 107, 1695 

Smith, B. J., Struck, C, Kenney, J. D. P., & Jogee, S. 

1999, AJ, 117, 1237 
Solomon, P. M., & Barrett, J. W. 1991, Dynamics of 

Galaxies and Their Molecular Cloud Distributions, 146, 

235 

Yoshida, M., Taniguchi, Y., & Murayama, T. 1999, AJ, 
117, 1158 

Young, J. S., Xie, S., Tacconi, L., et al. 1995, ApJS, 98, 
219 

Zhang, J. S., Henkel, C, Kadler, M., et al. 2006, A&A, 
450, 933 



Hunt et al.: NUGA: IX. NGC2782 



17 



CO(l-O) 




Fig. 1. ^^CO(l-O) velocity-channel maps observed with the PdBI in the nucleus of NGC 2782 with a spatial resolution 
of 2.1"xl.5" at PA=38° (beam is plotted as a filled ellipse in the bottom left corner of each panel). We show a field 
of view of 50", i.e. ^^1.2 times the diameter of the primary beam at 115 GHz. The phase tracking center is indicated 
by a cross at aj2ooo=09''14'"05.08^ and (5j2ooo=40°06'49.4". Velocity-channels are displayed from v=-240 km s~^ to 
v=240 km s~^ in steps of 10 km s~^. Velocities are in LSR scale and refer to v=2562 km s~^. Contour levels are -3a, 
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C0(2-l) 




Fig. 2. Same as Fig.[T]but for the 2-1 line of ^^CO. Spatial resolution reaches 0.7" xO. 6" at PA=20° (beam is plotted 
as a filled ellipse in the bottom left corner of each panel). We show a field of view of 18", i.e. ~0.9 times the diameter 
of the the primary beam at 230 GHz. Velocity-channels are displayed from v=-240 km to v=240 km s~^ in steps 
of 10 km s~^, with same reference as used in Fig. [1] Contour levels are -3(t, 3cr, 7a, 10a, 14a and 18a, where the rms 
(7=5.2 mjy beam^^. 



Hunt et al.: NUGA: IX. NGC 2782 



19 



o 

CD 
CO 
O 



<1 



-50 



1 1 1 1 


J 1 i 1 i 
•■ ♦ 


1 1 1 I 


1 1 1 1 1 
« 


1 1 1 1 1 1 


- 


■ 

« 

I 






* 


* 

- 




V 
■ 


m 


- 










- 

■ 












' 1 »■ 

1 . . . 


4 

, 1 , ,' . 


% 

1 .. . 


> . 

, , 1 , , 




150 


100 


50 





-50 






Aa (arcsec) 







100 

Aa (arcsec) 



0.5 



Fig. 8. Left panel: a) IRAC 3.6 /im image of NGC 2782. The stellar sheets to the east and west, and the ripples about 
~40" to the east are clearly visible. Right panel: b) IRAC 3.6 /im residuals from the bulge-to-disk decomposition 
described in Sect. 14.11 with contours superimposed. 
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Fig. 9. Left panel: a) shows CO(l-O) total intensity contours overlaid on the gray-scale residuals (in MJy/sr) of 
the bulge-disk decomposition performed on the 3.6 /xm IRAC image. Right panel: b) shows the analogous overlay of 
C0(2-l). The phase-tracking center is shown by a cross. 
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Fig. 12. Left panel: a) ^^CO(l-O) total intensity contoured over a false-color representation the IRAC 3.6 /im image. 
Unlike Fig. [71 both the IRAC and the CO images have been deprojected onto the galaxy plane and rotated (by 15°) so 
that the major axis is horizontal; units on [x, y) axes (AXg/AYq) correspond to arcsec offsets along the major/minor 
axes with respect to the AGN. The outer oval and stellar bar are marked by solid lines (in the deprojected rotated 
coordinate system). The image ficld-of-view is '^2.2 kpc in diameter. The (/)i-angles are measured from the +X axis in 
the counter-clockwise direction. Right panel: b) The same for ^^C0(2-l). 
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Fig. 13. Left panel: a) The different normalised Fourier components Qm (defined in Eq. [3l as the maximum tangential 
force, normalised to the radial force), for the potential of NGC 2782, derived from the NIR image. The top solid line is 
the total tangential force Qt, the next solid line is Q2, the dashed line Qi, the dotted line Qa, and the dotted-dashed 
one Qi. Right panel: b) The corresponding phases 9m in radians, defined in Eq. [U with the lines as in the left panel. 
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Fig. 14. Left panel: a) ^^CO(l-O) total intensity contoured over a false-color representation of the torque maps r 
derived as described in the text. As in Fig. I12[ this is the deprojected representation, with units on {x,y) axes 
(AXg/AYg) corresponding to arcsec offsets along the major/minor axes with respect to the AGN. The torque exerted 
on the spiral arm to the south is clearly associated with the outer oval. Right panel: b) The same for ^^C0(2~l). 
Both panels show the typical butterfly diagram due to an m = 2 perturbation, and in b) we have marked the four 
quadrants with their positive (I, III) and negative torques (II, IV) which drive the gas outward and inward, respectively. 
As described in the text, the sign of the torques is determined by the rotation direction of the gas, namely clockwise 
in NGC 2782. 




Fig. 17. Left panel: False-color plot of the stellar component distribution, in Run A (the intensity scale is linear 
and axes are in kpc). Middle: False-color plot of the corresponding gaseous component (the intensity scale is loga- 
rithmic). Right: Same as the middle panel, but with an expanded spatial scale. It can be seen that the gas which was 
predominantly in the more external ring (ILR of the primary bar) at T = 690 Myr, is falling progressively inward, and 
is primarily inside the ILR at T = 750 Myr. The gas infalling along the "dust lanes" is aligned along the nuclear bar, 
where the CO is observed in NGC2782. 
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Fig. 18. Left: Intensity of the m = 2 and rn = 4 Fourier components of the potential in Run A, during the second 
bar decoupUng. Right: Power spectrum of the m = 2 Fourier component, giving the pattern speed as a function of 
radius, during a period of 100 Myr corresponding to the coexistence of the two bars. Since the pattern speed of the 
nuclear bar varies with time during this period, its Ob signature is somewhat widened. 



